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Abstract: Recently a powerful duality between color and kinematics has been proposed for 
integrands of scattering amplitudes in quite general gauge theories. In this paper the duality 
proposal is extended to the more general class of gauge theory observables formed by form 
factors. After a discussion of the general setup the existence of the duality is verified in two 
and three loop examples in four dimensional maximally supersymmetric Yang-Mills theory 
which involve the stress energy tensor multiplet. In these cases the duality reproduces known 
results in a particularly transparent and uniform way. As a non-trivial application we obtain 
a very simple form of the integrand of the four-loop two-point (Sudakov) form factor which 
passes a large set of unitarity cut checks. 
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1 Introduction 

The ability to calculate is of fundamental importance in any would-be theory of nature. 
In high energy physics this manifests itself in the quest to calculate to higher orders in 
perturbation theory, including more and more loops and legs. The basic quantities to be 
calculated in a gauge theory are in general a mixture of plane-wave external states and 
gauge invariant operators. If only on-shell states are involved these observables are scattering 
amplitudes, while if only gauge invariant operators are involved these are pure correlation 
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functions. Admixtures of observables build of operators and on-shell states are called form 
factors. 

The main motivation of this paper is the observation that in recent years many new in- 
sights have been developed for the calculation of scattering amplitudes in gauge theories (see 
e.g. [1] and references therein), while similar developments for form factors and correlation 
functions have been minor in comparison. This paper is part of an effort to redress this imbal- 
ance. Recent developments have shown that maximally supersymmetric Yang-Mills theory is 
a good testing ground for new ideas which then should trickle down to non-super symmetric 
gauge theory. This line of thought will also be followed here. 

A particular insight that has recently been developed is a duality between color and 
kinematics that has been proposed for integrands of scattering amplitudes [2] [3]. This duality 
will be extended here to form factors. Color- kinematic duality in its current form states that 
at tree and loop level the integrand of a scattering amplitude can be arranged in such a way 
that whenever the color factors of an integrand obey a Jacobi relation, the kinematic parts 
of these diagrams can also be chosen such that they obey the same identity. In a second step 
which will not be of much importance for the purposes of this article, this duality can also 
be used to construct the integrand of a scattering amplitude in a gravity theory, whose field 
content is simply a tensor product of two gauge theory copies. The status of color-kinematic 
duality will be reviewed below. 

There are two distinct ways the ideas of color-kinematic duality could apply to more 
general observables which include gauge invariant operators. The first is simply as a direct 
consequence of color-kinematic duality for scattering amplitudes. On a unitarity cut and at 
tree level poles of quite general form factors the duality should apply to the amplitude parts 
of the residues and discontinuities, see for instance [4]. In this paper we conjecture that the 
duality can be used directly for general form-factors and in particular without taking any 
cuts. We will show in examples that the duality can be used as a powerful tool to produce 
consistent ansatz for high loop form factors. In particular, our results unify and greatly 
simplify the results in the recent literature up to three loops. 

As a new non-trivial application, we compute for the first time the four-loop integrand of 
the Sudakov form factor in AA=4 super Yang-Mills theory (SYM) in a form which passes all 
(highly non-trivial) unitarity checks performed. The existence of such a form for the four-loop 
form factor is a strong evidence that the color-kinematic duality holds in general. A distinct 
advantage of using color-kinematic duality is that the traditionally hard to obtain non-planar 
integrals in the problem are fixed by the coefficients of planar integrals. The evaluation of 
the found integrals is in progress [5]. 

This article is structured as follows. In section 2 the necessary background for this paper 
will be reviewed. This is followed by section 3 where our general calculational strategy is 
explained. This strategy is applied to a few examples in section 4, and in particular the four- 
loop form factor in section 5. The main presentation ends with a discussion and conclusion 
section. In the first appendix we give a representative selection of four-loop Jacobi relations. 
The other appendix briefly discusses color structures which appear in two-point form factors 
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and correlation functions with two adjoint operators through eight loops. 



2 Review 

2.1 Color-kinematic duality for scattering amplitudes 

Color-kinematic duality is the collective name for what is basically a pair of two partially 
proven conjectures [2], [3]. The starting point is to write the integrand for a scattering 
amplitude in D dimensions at I loops in terms of a sum over trivalent graphs only, 



where Cj is the natural color factor associated to trivalent graph i, Si is the symmetry factor 
associated to the trivalent graph i and Sj is the product of all internal propagators of the 
graph. This form of the integrand simply amounts to a rewriting of Feynman graph-based 
perturbation theory by parcelling out the four- vertex over three vertices. This determines 
some set of numerators rtj which is not unique. The non-trivial claim of color-kinematic 
duality is that there exists a set of numerators such that whenever color factors obey a Jacobi 
identity, the numerators do too: 



A set of numerators which obeys all Jacobi identities is called color-dual. Color-dual numer- 
ators have been constructed explicitly at tree level for all multiplicities through a variety of 
methods [6-8], and at loop level in specific examples in maximally supersymmetric Yang-Mills 
theory for five points up to two loops [9] and for four points up to four loops [10]. Moreover, 
in pure Yang-Mills theory numerators are known at two loops for the four-point helicity equal 
amplitude [3] and at one loop for the helicity equal and one-unequal cases [11]. Of course, 
the tree level results almost trivially will extend to any generalized cut of a loop amplitude 
which involves tree amplitudes. 

The second conjecture usually taken to be part of color-kinematic duality is that if a set 
of color-dual numerators exists in a gauge theory, then an amplitude in a gravitational theory 
may be constructed as 



where n and h are color-dual numerators for two in general distinct gauge theories. This 
conjecture has been proven at tree level, assuming that local numerators exist [12]. The field 
content of the gravity theory is the direct product of the field content of the gauge theories. 
In this way 'squaring' J\f = 4 super Yang-Mills gives the maximal (ungauged) N = 8 super- 
gravity theory. Similarly, squaring pure Yang-Mills theory gives Einstein gravity coupled to 
a dilaton and a two-form. The latter theory is sometimes referred to as M = supergravity. 





V{i, j, k] Ci + Cj + Cfe = =>• rii + rij + rife = 



(2.2) 




(2.3) 
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The color-kinematic duality relation of gauge theories to gravitational theories will not play 
a central role in this article. 

As an illustration of the general idea, consider the four-point tree amplitude in gauge 
theory which can be written in terms of three trivalent diagrams as 




Figure 1. Trivalent diagrams for four-point tree amplitudes. 



The color factors denoted as q are given by the product of f abc associated to each trivalent 
vertex 

Cs — Jaia2b j6<23a4 ^ _ ja2asb jba^ai ^ _ !ja\a-zb ^baia^ g\ 

where 

f abc = iV2f abc = Tx{[T a ,T b ]T c ) . (2.6) 
They satisfy the Jacobi identity 

c s = c t + c u . (2.7) 
By color-kinematic duality one can also have 1 

n s = n t + n u . (2.8) 

The solution to the Jacobi conditions may be given explicitly as 

n s = a, th = — ^ ^£ 0) (1, 2, 3, 4) -a-, n u = t Af ] {1,2,3,4) - , (2.9) 

This can be shown to reproduce the full color-dressed amplitude. The solution depends on 
an arbitrary parameter a. This freedom in specifying numerators is the so-called generalized 
gauge transformation. This gauge freedom is the freedom to choose another set of kinematic 
numerators 

< = Tn + Ai, (2.10) 



Note that we have necessarily made some choice of signs for the color factors. One can choose other 
conventions, and the relation between rii should change correspondingly. The final result is given in terms of 
the product a m, and therefore is independent of this choice. 
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Figure 2. Color-kinematics duality at loop level. 



without changing the tree amplitudes as long as 

^^i = o (2.11) 
r 

Moreover, the A should satisfy Jacobi identities. 

It should be obvious that through unitarity cuts tree-level color-kinematics duality pro- 
vides constraints to the loop integrand. This occurs whenever a tree amplitude is isolated by 
the cuts. For example as shown in Figure 2, if we cut four propagators £±, ..,£4, we have 2 



A® 



(£f,..,^)-cut 



JdUPSih, .., £ 4 ) 4+4 4) 4 0) (4,^2,4,4) 



c s n s c t n t c u n u 



(h+l2? (£i+£ 3 ) 2 (h+£ 3 ) 2 _ 

(2.12) 



Cj,sNj^ s CjjNjj Cj,u^j,u 

h + h) 2 (h + Q 2 (h + hf 



where in the second line the BCJ representation of four-point tree amplitudes of equation 
(2.4) was used and ^,4' has been expressed in terms of some loop integrals. The third line 
is simply given by defining 

Cj >a = cj c s , N j)S = rij n s . (2.13) 
Because of n s = n t + n u , we have 

which may be thought of as the dual relation of that of color factors, Cj jS = Cjj + Cj tU . 
Beyond amplitudes it was noted in [4] that a similar observation should hold for unitarity 
cuts 3 of general form factors, whenever a scattering amplitude is isolated on the cuts. 

Beyond the constraints on loop integrands from unitarity cuts it was conjectured that 
there exists a representation where the above duality relations are true off-shell, i.e. without 



2 Sums over all states to appear in the cuts are understood but not shown explicitly in the formula. 



3 This includes the case of tree level poles. 
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doing cuts to impose on-shell conditions [3]. More explicitly, given a representation that 

i i 3 

there is an off-shell color-kinematics duality 

C l = C j +C k => Ni = Nj + N k . (2.16) 

This is a highly non-trivial claim. This has been checked on a case-by-case basis, but no 
general proof exists. Some supporting evidence was recently uncovered in [13]. It may be 
taken as an evidence of the existence of a Lagrangian formulation of gauge theory which has 
explicitly color-dual structure. This much more powerful version of color-kinematic duality 
at loop level will be generalized to a class of form factors in this article. 

2.2 Form factors in M=4 supersymmetric Yang-Mills theory 

As recalled in the introduction, form factors are the most generic class of observables in any 
Yang-Mills theory. They are mixtures of i gauge invariant operators O(x) and j on-shell 
states s, 

F(0 1 ...O i ,s 1 ,...s j ) (2.17) 

In this article only the simplest form factors with one insertion of a gauge invariant operator 
will be considered. This will be referred to as a j-point form factor. The simplest example is 
the two-point form factor which is also called the Sudakov form factor. In this article these 
form factors are considered in the context of maximal N = 4 super Yang-Mills theory. The 
operator inserted will belong to the half-BPS stress-energy multiplet: the multiplet containing 
the (on-shell) stress energy-tensor as well as the self-dual and anti-self-dual field strengths. 
Another member of the multiplet which is usually useful for practical calculation is 

O = Tr(</> 3 V 34 ) (2.18) 

where the field (j) IJ with anti-symmetric SO (6) indices parameterizes the six scalars of the 
M = 4 super Yang-Mills theory. 

The object of study in this article is the form factor with the operator Fourier transformed 
to momentum space, i.e. 

F(l,...,n;q) = J d 4 x e~ iqx (1 • • • n\O(x)\0) , (2.19) 

The momentum of the inserted operator in general does not square to zero, q 2 ^ 0. In this 
sense it is 'off-shell'. This does not mean it does not satisfy field equations: the conservation 
of the stress-energy tensor in field theory for instance is typically only guaranteed up to field 
equations. For the form factor these are satisfied for all inserted operators in the problem. 

Compared to progress on amplitude computation, form factors are much less studied and 
have only attracted renewed attention very recently, see for example [14-24]. The current list 
of achievements at weak coupling beyond one loop is 
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• 2pt @2 loops: [14], 

• 2pt @3 loops: [22] 

• 3pt @2 loops: [23] 

Moreover, the 3- loop result was only obtained after the much harder QCD computation [25]. 
Form factor at strong coupling (in 2D kinematics) was considered in [15, 16]. In the following 
mainly two and three points form factors will be considered, unless indicated otherwise. These 
have both special kinematics as well as special color properties. The two-point form factor has 
to be proportional to the color delta function and can only be a function of the momentum 
squared of the operator leg, q 2 as well as, in SU(N C ) gauge theory, the number of colors N c 
and the gauge coupling g ym . Moreover, the two-point form factor has to be proportional to 
the tree form factor by on-shell super symmetry, 



In particular, the function / does not depend on the inserted operator. As a function of N c 
and g ym it can be rewritten in terms of powers of the 't Hooft coupling, 



and, in general at a finite loop order as a finite taylor series in N~ . Non-trivial terms in this 
series are called non-planar corrections. Actually, as far as the color structure is concerned, 
the two-point form factor is the same as the two-point correlator which is computed to extract 
the (3 function at high loop order [26] . From this it follows that the first non-planar correction 
to the two-point form factor occurs at four loops, regardless of gauge theory under study. This 
point is worked out more fully in appendix B. 

For general gauge group the number of colors N c in an SU{N C ) calculation is replaced 
by appropriate Casimir operators. The three-point form factor can have more interesting 
dynamics and more complicated color factors. In particular, in addition to the structure 
constant for the gauge group the completely symmetric structure constant d abc might make 
an appearance. On-shell supersymmetry still restricts the form factor to be proportional to 
the tree. Since the latter is known to be absent at two loops [23] in M = 4 SYM, the first 
appearance of this can be at three loops. Based on the details of the computation of the beta 
function at three loops in [27-29] it seems likely it indeed does appear at this loop order. 
For more points more general color structures and operator-dependent pre-factors can and 
do appear. 

3 General formalism 

In this section first color-kinematic duality will be generalized to form factors. This is first 
done at tree level for clarity and then extended to loop level. In the second half of the section 
we outline how this conjecture can be used in general to compute the integrand of the form 
factor. This technique is then applied in the following sections in specific examples. 



F«(l, 2; q) = 5 ab F^(l, 2, q)f{q\ N c , g ym ) 



(2.20) 



(2.21) 
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3.1 Color-kinematic duality for form factors 



The color structure of form factors is different from that of amplitudes. To get our bearings a 
few explicit examples will be studied first. The simplest example is the two-point case which 
always has the trivial color factor 5 ab where a and b are the color indices of the gluons. In 
particular, there is no Jacobi identity This also holds for the three-point form factor at tree 
level as it can be written as the sum of three trivalent diagrams as 



(0) 



JF. 



where 



ja^a-i <5(8)(^ ir?1 + x 2m + A 3 t? 3 ) _ f aiCl2b 5 asb n l2 + f a2asb S aib n 23 + f a3aib 5 a2b n 31 



(12) (23) (31) 



Sl2 



S23 



S31 



71; 



-<5 (8) (Ai7]i + A 2 ?72 + A3r ? 3) 



(3.1) 
(3.2) 



v 3(12)(23)(31) 

The fermionic delta function is included, see e.g. [17, 18]. Hence this form factor is propor- 
tional to a single structure constant. 

The four-point case becomes more non-trivial. There are 15 different trivalent diagrams: 
there are five ways the gauge invariant operator can be inserted in the trivalent graphs of 
a four particle amplitude. Some of the graphs are illustrated in figure 3. On poles of the 
Feynman graphs where a tree level amplitude is isolated color kinematic duality should hold. 
Consider for instance the sum of the three diagrams shown in Figure 3(a) 



0(4) 



£0,4 c 

S123 



ja\a,2b jbazc 



fa 
12 | J 



0,20.3b jbcai ^,(4) faiasb fba2C ^Vy 



n 



fa 

23 1 J 



f" 



n 



13 



S12 



S23 



S13 



(3.3) 






(a) 



(b) 



Figure 3. Diagrams for four-point tree form factor. 



At the kinematic pole S123 = the residue given by the expression in the bracket becomes a 
four-point amplitude. Color-kinematic duality for the tree amplitude then imposes 

4 , L,^„ = (4 , + 4t , )l„^»- ("J 

This is the lesser sense in which color kinematic duality should hold for form factors. The 
conjecture we make in this paper is that the duality also holds before taking any poles, in a 
sense to be made precise below. 
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Before writing more general formulas, let us illustrate the idea in the case of the tree 
level four-point form factor. In this case the form factor can be written over two classes of 
diagrams: those with the operator inserted on an external leg (class D, see equation (3.3)) 
or the class with the operator inserted on the internal leg (E). 



t[ 0) = J2d^+E (3.5) 



-(0) _ 

i=l 

Explicitly, E follows from the diagrams such as shown in Figure 3(b) 



j&c/W^W jbc Ja 2 a 3 b Ja^c gbc f ai a 3 b fa 2 a 4 c 

£> = 1 1 , (,3.bJ 

S12 S34 S23 S41 S13 S 2 4 

We conjecture that numerators can always be found such that they mirror the Jacobi identities 
of the color structures they multiply, i.e. 

41 = 41+41, (3.7) 

as well as 

n (12,34) = n (23,41) + n (13,24) > (3-8) 

This is quite plausible as a simple counting argument shows. At tree level, there is 
a minimal basis for color factors modulo all Jacobi relations as explained in [30] (see also 
[31]). For four particle color factors ( "colormatics" ) this minimal solution has two elements. 
Implementing this for both color as well as kinematic factors in equation (3.5) gives two 
equations in terms of 10 unknowns, i.e. 2 for each D and 2 for the E type factors. This 
set of equations has many solutions. In fact, it is straightforward to check that the resulting 
equations have solutions with additional symmetries. In particular, one would expect that 
the different D type diagrams are related by simple interchange symmetries. 

One solution is to set all E type numerators to zero and 

"12 = \ 0i2(s23 + 524)^(1234) + Si2 s 24j F(1324)) (3.9) 

"23 = \ ( S 24(S23 + 812)^(1324) + 512*24^(1234)) (3.10) 

where 

Sij = (Pi +Pj) 2 with J *' (3-11) 

Pi = p 4 + q , 

and F(1234) as well as F(1324) are color ordered form factors with the indicated ordering. 
The other nW are simply defined by interchanging 4 with i in the above formula. This solution 
was obtained by explicitly solving the system of equations. 

In particular, in this solution it is seen that no constraints on the form factors are nec- 
essary: there are no BCJ-type [2] relations for the four-point tree level form factor. This 
is to be contrasted with the four-point amplitude computation where a similar computation 
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shows that the duality imposes additional relations on the amplitudes. The color kinematic 
duality is expected to only hold for amplitude subgraphs separately. The sets of color-dual 
numerators do not in general sum to anything physical separately. 

The above observation for four-point form factors can be generalized to form factors at 
both tree and loop level. At tree level, one should decompose the tree level form factor in 
terms of the place the gauge invariant operator is inserted, 



4 0) = E^ } '° ( 3 - 12 ) 



Then for each J-^' separately we conjecture color-dual numerators can be found. There is 
a considerable number of freedom to choose the numerators as a simple counting argument 
of equations and unknowns will yield. 

Color-kinematic duality for loop form factors 

At loop level a similar decomposition may be imposed. Then for every one introduces 

numerators and color factors as in the amplitude integrand case, 

4^ = £/n<^r (3.i3) 

r\ J 3=1 d * 4 

Moreover, one can impose numerator symmetries as in the D-type contributions as above: if 
two diagrams are related by a symmetry, then their numerators are related. At tree level, by 
an extension of the counting argument above it is quite plausible that color-dual numerators 
always exist. At loop level it will be shown below that the conjecture gives strong constraints 
for the integrand of form factors. 

This conjecture at loop level will be explored below in several examples in the next section. 
In the remainder of this section a complete overview is given over the calculational setup. 
The basic outline of this mirrors to quite some extent the setup for scattering amplitudes up 
to four loops as described in [10]. 

3.2 Calculational strategy 

The full integrand of a L-loop form factor can be given as 

4 L ) =K n Y,Y,T. CiIi > (3 - 14) 



On I 



where we sum over a set of basis integrals 4 

L 



4 The coupling constant g ym will usually be suppressed in our results. 
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which are associated to L-loop diagrams with only trivalent vertices. The classification of 
trivalent graphs will be discussed in the following subsection. We sum over all permutations 
a n of external on-shell legs. K n is some uniform kinematic factor. For two-point and three- 
point cases that we will consider 

K 2 = sj 2 4°\ K 3 = 4 0) . (3.16) 

The Cj's are color factors. For a given trivalent graph, each cubic vertex can be dressed 
with a factor f abc . In the form factor there is also a special vertex, the one connected to 
the off-shell leg q that is dressed with the factor 5 a b = Tv(T a T b ). The whole color factor d 
is then given by taking the product of all f abc, s and the single 5 a b- The Si's are symmetry 
factors. They account for the over-counting both from the sum over the permutations of 
external on-shell legs and the internal symmetries of the graph. In practice, Si for a given 
graph can be computed by randomly shuffling the list of edges and then counting the number 
of inequivalent isomorphic maps. 

The main challenge is to determine the numerator iVj for each basis integral. In principle, 
one can apply Feynman rules and sum over all Feynman diagrams. However, this becomes 
very complicated for higher loops, in particular for AA=4 SYM with its large number of fields. 
An alternative is to construct the amplitude based on traditional unitarity methods. 

Here the conjectured color-kinematics duality for the loop integrand is used to construct 
a consistent ansatz which is then verified through unitarity cuts. For a given trivalent graph, 
one can apply the Jacobi relations to relate it to two other graphs, such as those shown in 
Figure 2 

Ci = Cj+C k =► Ni = Nj + N k . (3.17) 

Note that no cuts for loop propagators are imposed here. For each propagator (except the 
two connected to the off-shell leg q), there is one such equation. By considering all graphs, 
one can obtain a set of equations for the numerators. This gives strong constraints on the 
numerators. The main job is to find a solution which satisfies the full set of equations, while 
also consistent with the unitarity cuts. Below the various substeps in the calculation are 
discussed in yet more detail. 

3.2.1 Generating topologies 

The first step in applying color-kinematics duality is to generate all trivalent graphs at a given 
loop order. There are mature techniques to generate inequivalent topologies; we have opted 
for [32] because of its simplicity and ability to handle higher loop orders. Manipulations of 
graphs have generically been performed with Mathematica. For the problem at hand we need 
to make some particular choices. 

1. The diagrams are trivalent, one-particle irreducible graphs, without bubble sub- graphs. 
For the form factor graph the vertex which is connected to the off-shell leg q (corre- 
sponding to the operator) is special. 
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2. Since we mainly focus on jV=4 SYM, we further exclude most diagrams which contain 
triangle subgraphs, due to the no-triangle property of AA=4 SYM [33]. For form factors, 
a single triangle is allowed if it is connected to the off-shell leg q. The triangle subgraphs 
we exclude are the one-loop and two-loop non-planar triangles 5 . This is enough up to 
four loops. For five and higher loops, there can be three-loop or higher-loop triangle 
subgraphs, which may be excluded and can reduce further the number of topologies. 
(This further subtraction is not considered in the following table.) 

3. The diagrams usually have the maximal number of loop propagators. In the construction 
with color-kinematics duality, one also needs to consider the graphs with less loop 
propagators but with more complicated trivalent tree legs, for example, the "snail-like" 
diagrams [10]. Such diagrams can be easily obtained by replacing a single leg with multi- 
leg trivalent tree graphs in lower-point cases. One could also generate these graphs 
by applying color-kinematics relations to those graphs containing maximal numbers 
of propagators. (In the following table only graphs with maximal numbers of loop 
propagators are counted.) 

Based on the above selection rules, we obtain the number of different topologies for two 
and three points form factors and four and five-point amplitudes up to six loops, as shown in 
Table 1. 

Table 1 . The number of topologies of trivalent graphs allowed by the selection rules. 



# of loops 


1 


2 


3 


4 


5 


6 


# of topologies for A± 


1 


2 


9 


52 


446 


4891 


# of topologies for F2 


1 


2 


6 


34 


273 


2718 


# of topologies for A$ 


1 


3 


19 


155 


1684 


22225 


# of topologies for F3 


1 


4 


22 


171 


1695 


20046 



From the table it follows that for higher loop two-point form factors the number of graphs 
grows about as fast as for the four-point amplitude, but is typically smaller than the four 
point amplitude by a factor of 2. 



5 For two-loop non-planar sub-triangles this concerns graphs like 
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3.2.2 Solving color-kinematics dualities 

After generating the trivalent graphs, the constraints imposed by color-kinematic duality can 
be studied. 



1. Given a trivalent graph, for each internal propagators (except the two connected to the 
off-shell leg q), one can apply the Jacobi identity to generate two other graphs. These 
two graphs can both be part of the basis set. Often, one of the two contains a triangle 
subgraph. In this case the corresponding numerator is set to zero. By considering all 
propagators of all graphs, one generates a set of equations for the numerator functions 
of the trivalent graphs in the basis. 

2. Typically the numerators of a small set of graphs can be identified from which all other 
numerators can be obtained by the set of equations derived from the duality. The 
corresponding graphs for these numerators are called master integrals. These can be 
identified by inspecting the full set of equations: there is usually no unique choice. 

3. One can construct the most general ansatz for the master integrals, under several simple 
constraints derived from physical expectations: 

• One would expect the numerators to be local. In particular they should be poly- 
nomial functions of momenta (both loop and external). 

• The full expression should not break the excellent UV properties of M = 4 SYM. 
In practice this implies that for any n-point one-loop subgraph, the numerators 
contain no more than n — 4 powers of the loop momentum for that loop. The 
exception is that if the one-loop subgraph is a one- loop form factor, one allows the 
maximal power to be n — 3. 

• The numerator should preserve the symmetry of the corresponding graph. 

• The numerator should be consistent with the maximal cut (where all internal legs 
are cut). One can apply the simple "rung-rule" in such cuts [34]. 

By applying these constraints one obtains an ansatz for the master integrals with a 
number of free parameters. 

4. Given the ansatz for master integrals, one can obtain the numerators of all other inte- 
grals in the basis by using a small set of Jacobi relations. We then impose the same 
constraints considered above for all integrals. This typically fixes a large number of 
parameters. It is then necessary to check if the obtained numerators satisfy the full 
set of color-kinematics equations. If all Jacobi relations are satisfied, we have found a 
color-dual solution. 

5. It is not guaranteed that the solution obtained is physical and in general free parameters 
are left at this stage. To both fix these parameters and to make sure that the solution of 
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the integrand is indeed physical, one should check it against a complete set of (preferably 
D-dimensional) unitarity cuts [33, 35]. In the next subsection, we discuss how to do 
these unitarity checks in a systematic way. 

3.2.3 Unitarity checks 

On a given unitarity cut one can compute the cut integrand both from the ansatz as well as 
the product of lower loop amplitudes and form factors, summed over the complete spectrum 
of the theory under study. For the solution obtained from the procedure outlined above the 
two expressions must be consistent with each other i.e 

k » e e j. c * * L = / n ^ Ftrcc n -< rce . (3.i8) 

On i 1 J i a 

In practice it is convenient to consider color-stripped cuts involving tree amplitudes. The tree 
inputs are color order amplitudes or form factors, and the cut integrand on the LHS can be 
also constructed without considering color factors. 

Although the idea is straightforward, in practice there are two technical problems. One is 
to compute the product of tree amplitudes, where the main challenge is to do the summation 
over the on-shell states appearing on cut legs. The other is to extract the cut integrand from 
the ansatz, in particular when there is a large number of basis integrals. 

Tree input 

We can start from supersymmetric expressions for tree amplitudes and form factors, and take 
the product of them by integrating out the fermionic variable 7]g. for each cut leg £{. For 
example, in the four-loop computation, we need to consider a quintuple cut which needs the 
following tree input 

r 5 r 

J II d ^k [^5 MHV Hi, -h, ~h, ~k, ~h) Af KV ( Pl ,p 2 , h,k, h,h, h) (3.19) 



i=l 



jrNMHV^_^ MHV^^p^ ^ j 4) ^ ^ 

jrN mhv^_^ A^ Mm/ (pi,p2,h, h, h, h, h) 

jrmax-non-MHV^_^ ^MHV^^^ J g) / 4j / 3) \^ ^) 

The main complexity is the fermionic integration as it involves high MHV degrees. This has 
been discussed for example in [36, 37]. In practice we find it convenient to apply the MHV 
rules as in [36]. The supersymmetric form factor was developed in [19, 20]. 

Cut integrand from the ansatz 

For the higher-loop integrand given by a large number of integrals, there would be hundreds 
of cut diagrams contributing to some unitarity cuts. For example, for the above quintuple 
cut of the four-loop form factor, there are 192 cut diagrams to sum over. Obviously, one 
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has to develop a systematic method to automatize the computation. In principle, one could 
start from the basis graphs, consider various possible cuts of them, and then collect the 
contributions. In practice, we find the following way more convenient. A similar idea has 
been used in [38]. 

1. Rather than starting from the basis and then doing cuts, we "inverse" the procedure. 
We first construct the trivalent tree diagrams of amplitudes and form factors. We take 
them as the cut tree diagrams in the RHS of (3.18). 

2. Given the trivalent tree graphs, the next step is to construct loop diagrams by sewing 
trivalent tree diagrams in all possible ways. Since the tree diagrams are trivalent, the 
loop diagrams obtained are also trivalent. Not all such loop diagrams correspond to 
basis integrals. If a loop diagram is isomorphic to one of the basis, then it corresponds 
to a cut contribution of that basis. If it is not isomorphic to any basis, we can simply 
neglect it. In this way, by sewing all possible trivalent diagrams, we are guaranteed to 
obtain all possible cut contributions from all basis integrals, without double counting. 

3. The cut integrand for each cut diagram is given by the product of the corresponding 
numerator and tree propagators. In particular, the numerator Ni for a given cut diagram 
should be re-expressed in terms of external and cut momenta. As we mentioned that it 
is convenient to consider color stripped cuts, it is not necessary to consider color factors 
here 6 . The full cut integrand for a given cut is then given by summing over all cut 
diagrams 

K n Ni (tree propagators) , (3.20) 

cut graphs; 

which should be compared with the integrand obtained from the product of tree ampli- 
tudes. 

The above procedure has been applied to a few known results up to three loops and all 
found consistency, which shows that the method is robust. We then apply it to the four- loop 
form factor, such as the above quintuple cut of four-loop form factors, which also provides 
consistent results. 

4 Examples 

In this section the general technology outlined in the previous section is applied to various 
form factors in maximally supersymmetric Yang-Mills theory above one loop 7 which have 

6 We add a comment on the non-planar cut, which can be obtained by permuting the external legs of tree 
amplitudes. For these cases it is crucial to obtain the correct sign factor for each cut diagram that is related to 
the order of tree legs. One way to obtain these sign factors is to assign a color factor for each trivalent vertex 
of tree graphs, and compare the "glued" color factor with the color factor of the basis. 

7 It is easier to work out the one-loop two and three points form factors following the higher-loop examples 
considered here. 
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appeared in the literature. Although the results are known, our techniques illustrate the 
existence and the power of the duality. 

4.1 Two-point two-loop form factor 

As a warm-up exercise, we consider first the two-loop two-point form factor. This result has 
been computed by Feynman graph methods in [14]. 




(a) (b) 



Figure 4. The integrals for the two-point two-loop form factor. 



First, by equation (2.20) the two-point form factor in maximally supersymmetric Yang- 
Mills theory is trivially dependent of the inserted operator through a tree factor, as long as 
it's in the stress energy tensor multiplet. This tree factor is factored out, as will be done in 
every two-point calculation in this article. 

By the rules introduced in the previous section, there are only two trivalent graphs to 
consider as shown in Figure 4: a planar ladder and a non-planar ladder diagram. The Jacobi 
relations simply tell us that the numerators of both integrals are the same. 

By the power counting constraint explained in the previous section, we find that the 
numerator should be independent of loop momenta. Hence by the kinematics of the problem 
the numerator should be proportional to a power of s±2, which will be absorbed into the 
whole kinematic factor = sf 2 ^2 • The numerator is then a purely numerical constant. 
This numerical constant can be easily fixed by considering any (color stripped) unitarity cut, 
which turns out to be one. At the same time, this unitarity cut verifies that nothing has been 
missed in the construction. In this simple example it is not hard to explicitly compute and 
verify all possible unitarity cuts, verifying that the result is physical. 

The results including the color and symmetry factors are summarized in Table 3. The 
full form factor result can be obtained as 

b 1 

(T2 i=a 

= 5 aia * s\ 2 JT(°) (4 h + I 2 ) , (4.2) 

which reproduces exactly the known result [14]. Note that the color and symmetry factors 
are responsible for the numerical integer factors, which are 4 and 1 for planar and non-planar 
graphs respectively. 
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Table 2. The result for the two-point two-loop form factor. 



Basis 


Numerator factor 


Color factor 


Symmetry factor 


(a) 


1 


4iV c 2 (5 aia2 


2 


(b) 


1 


2 AT? 5 aia2 


4 



4.2 Two-point three-loop form factor 

As a more non-trivial example, the two-point form factor at three loops is calculated next by 
the procedure outlined above. This result has been computed by unitarity methods in [22] . 

First, by generating topologies we can find there are six trivalent diagrams, as shown in 
Figure 5. 8 . 




(d) (e) (/) 



Figure 5. The integrals for the two-point three-loop form factor. 



By applying the color-kinematic relation to this set of trivalent diagrams, a set of equa- 
tions can be obtained for the numerators. It turns out that one can choose the single integral 
(d) as the master integral. One can then make an ansatz for the numerator of this master 
integral by applying the following three constraints. 

1. From the power counting property, the numerator should depend only linearly on the 
loop momentum I and there should be no dependence on other loop momenta. A general 
ansatz is therefore given as (note that we have factorized a whole factor sf 2 ) 

Nt satz (Pi,P2j) = ail ■ pi + a 2 £ ■ P2 + a 3Pl • p 2 , (4.3) 

8 There is one bubble-like graph containing a two-point tree leg which turns out not contribute. For simplicity 
we do not include it here. In the four-loop construction, such graphs are shown as in Figure 10. 
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which contains three parameters ai,i = 1,2,3. 

2. The numerator should preserve the symmetry of the graph, which implies that it should 
be invariant under 

{pi,p 2 ,£} {P2,Pi,q~i}, (4.4) 

or more explicitly 

N d (pi,p 2 ,£)=N d (p 2 , Pl ,q-£). (4.5) 
Plugging in the ansatz (4.3), we obtain the relation 

a 2 = —a\ . (4-6) 

3. Finally, we consider the constraint of maximal cut. From the "rung-rule", we read the 
numerator {£ — p\) 2 . On the maximal cut we have 



[N d ( Pl ,p 2 ,£)-(£-pi) 2 ] =0. (4.7) 

maximal cut 

This fixes the remaining two parameters 

ai = -1, a 3 = -1 . (4.8) 

Therefore, by applying the above constraints we arrive at a unique solution for the master 
integral 

N2r** = (p 2 -p 1 ).£-p 1 .p 2 . (4.9) 

Given this solution, one can check that all Jacobi equations are satisfied. Other numer- 
ators can be obtained from the master integral by using the relations 

N a = N b = N c , N d = -N e = N f , (4.10) 
N b {pup 2 ) = -N e ( Pl ,p 2 ,£) - N e (p 2 , Pl ,£) , (4.11) 

where N x = N x (p±,p 2 ,£) if not specified. 

Now it is essential to check that the solution is indeed physical i.e. satisfies the unitarity 
cuts. A non-trivial quadruple-cut is given as in Figure 6. The product of trees is 

/ [J A*« [-F 4 MHV Hi, -l 2 , -l 3 , -l 4 ) A™ HY ( Pl , P2 , h, k, l 2 , h) (4.12) 
J i=i 

jtNmhv^_^^ Aq (pi,p 2 ,h, hih,h) 

^max-non-MHV^^ _^ _^ _^ ^MHV^^ ^ ^ ^ ^ 

From the basis integrals, we obtain the cut integrand as a sum of 29 cut diagrams. We have 
compared the two expressions numerically, and have found perfect agreement. 
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The results have been summarized including the color and symmetry factors in Table 3. 
The full form factor result can be obtained as 

^ 3) = ^ 0) ££^- ( 4 - 13 ) 

(72 i=a 

Note that the graph (f) has zero color factor and therefore does not contribute to the final 
result of the form factor. However, it is necessarily involved in solving the Jacobi relations. 



Table 3. The result for the two-point three-loop form factor. 



Basis 


Numerator factor 


Color factor 


Symmetry factor 


(a) 


s 2 
s 12 


8iV|<5 aia2 


2 


(b) 


S 2 
6 12 


4iV 6 3 (5 aia2 


4 


(c) 


S 2 
b 12 


4iV| <5° ia2 


4 


(d) 


(?2 -Pi) -i-pi -P2 


2iV| <5 ai ° 2 


2 


(e) 


~(P2 -pi)-t + Pl-P2 


2N?5 aia2 


1 


(f) 


(P2 - Pl)-£-pi-P2 





2 



The result we obtain by applying color-kinematic duality seems quite different from that 
in [22]. However, it is a simple check that the results are equivalent, by using the identities 
given in section 3 of [22] between different integrals. Our result is presented in a much simpler 
form which involves only trivalent graphs. The numerical integer factors which have no easy 
interpretation in [22] are also naturally explained here by the color and symmetry factors. 
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4.3 Three-point two-loop form factor 

As an example containing more external legs, we present the three-point two-loop form factor 
in a form satisfying all Jacobi equations. The three-point two- loop form factor was computed 
in [23]. 

There are six trivalent diagrams, as shown in Figure. 7. 




(d) (e) (/) 

Figure 7. The integrals for the three-point two-loop form factor. 



By inspecting the set of Jacobi relations, we can find that one needs at least two master 
integrals, which may be chosen as graph (a) and graph (d). The numerators of the other 
graphs can be obtained by the following relations 

N b (pi,p 2 ,P3) = N a (p 2 ,P3,pi) + N a (p3,P2,Pl), (4.14) 

N e (pi,P2,P3,Z) = N a (p 1 ,p 2 ,p 3 ) - N d (p 1 ,p 2 ,p 3 J), (4.15) 
N c = N b , N f = N d , (4.16) 

where N x = N x (pi,p2,P3, £) if not specified. 

Using the numerators in [23] for two master integrals, one can check that all the Jacobi 
relations are automatically satisfied. The numerator solution is given in Table 4. The color 
and symmetry factors are also given in Table 4. The full result can be constructed from the 
table as 

^ 2) = ^ 0) EEj«- ( 4 ^) 

1 J i 

It is a simple check that the result is equivalent to the one presented in [23]. We have also 
checked that the solution passes various triple-cuts following the procedure of section 3.2.3. 

5 Four-loop two-point form factor 

In this section the full four-loop integrand of the two-point form factor is constructed. To our 
knowledge this has never been computed before. 
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Table 4. The result for the three-point two-loop form factor. 



Basis 


Numerator factor 


Color factor 


Symmetry factor 


(a) 


Q 2 «12S23/2 


2 jy2 ja\a2a-i 


2 


(b) 


<? 2 S12(S13 + S23)/2 


1 N 2 f aiCL2a 3 


2 


(c) 


q 2 Sl2(si3 + S 23 )/2 


2 jy2 Ja 1 a2az 


4 


(d) 


Sl2(sis£-Pl ~ S 23 £-p 2 ) 


jy2 ja 1 a2a s , 


1 


(e) 


q 2 S\2S2z/2 - Si 2 (si 3 £ • Pi - S23^ " P2) 


jy2 ja 1 a2a s 


2 


(f) 


Sl2(si 3 £-pi - S 23 £-p 2 ) 





2 



5.1 Outline of the construction 

The construction is similar to the three-loop example. Here we outline the main steps. 

1. There are in total 37 trivalent diagrams to consider, which are as given in Figure 8—10. 
For completeness, we consider also three bubble-like graphs as shown in Figure 10. 

2. By applying the Jacobi relations for all possible propagators of all 37 integrals, one 
obtains a set of equations for the numerator functions. Part of the equations are given 
in appendix A. From the equations, one can identify two integrals as the master integrals: 
(13) and (21), which both have planar topology. 

3. We then construct an ansatz for the numerators of the master integrals. By power- 
counting, one can write down a most general ansatz as a polynomial of loop and external 
momenta for the numerator. For the master integral (21), the numerator depends at 
most quadraticly on £ 3 , so a general ansatz is 

N 2 i = ai(4 - pi) 2 + a 2 (l3 -P2) 2 + "3(^3 -Pi)(h -P2) 

+ [04(^3 • 4) + 05(^3 • Pi) + a 6 (l 3 ■ p 2 ) + a 7 (pi ■ p 2 )] (pi ■ p 2 ) . (5.1) 

For the master integral (13), the numerator depends linearly on £ 3 and £4, and a general 
ansatz is 

N 13 = /M4 • Pl)(h ■ Pi) + &(4 • P2)(h • P2) + • Pl)(h ■ P2) + /3 4 (4$ ■ P2){h ■ Pi) 

+ [&(4 • h) + fo(l3 ■ Pl) + fo{h ■ P2) + ■ Pi) + A>(4 ■ P2) + 010 (Pi • P2)} (Pi ■ P2). 
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(1) 



VI) 



(3) 



(-1) 




(17) (18) (19) (20) 



Figure 8. The basis (l)—(20) of two-point four-loop form factor. These diagrams have only planar 
contributions. Six of them have zero color factor. 

There are in total 17 parameters. The numerators of all other basis elements can be 
determined from these two via Jacobi relations. 

4. Then we can impose the constraints from the maximal cut 9 . There are five planar 
9 We apply the maximal cut and symmetry constraints for master integrals and other integrals together. 
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(21) (22) (23) (24) 




(31) (32) (33) (34) 



Figure 9. The basis (21)— (34) of two-point four-loop form factor. These 14 integrals have non-planar 
contributions. Five of them also have planar contributes. 




(35) (36) (37) 



Figure 10. The three bubble integrals which do not contribute to the final result. 



integrals for which one can easily read off the numerators directly from the rung rule: 



iVf 


— s 12 


3 


(5.2) 


N£ 


= S12 


(h -Pi) 2 , 


(5.3) 




= S12 


(h-hf, 


(5.4) 




= (h 


-Pif(h-Pi -P2) 2 , 


(5.5) 




= (4 


-Pi) 4 - 


(5.6) 
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In the maximal cut the ansatz must be consistent with these coefficients i.e. 

(jWi ^Vj ) | maxima] cu t ^ ' (5'^) 
These equations fix ten parameters. 

5. The numerators should satisfy the symmetry of the diagrams. Since the numerators of 
many graphs are related to each other by Jacobi relations, it is enough to consider the 
symmetries of the following graphs 

{(1), (6), (9), (12), (20), (21), (25), (26)} . (5.8) 

For example graph (26) has the symmetry: 

{pi,P2,4j^4,4s,M <-> {p2,pupi2 -4,4 + 4-4,4,4}- (5.9) 

These symmetry constraints can fix 4 parameters. Together with the above constraints 
from the maximal cut, one can fix in total 13 parameters, which can be explicitly given 
as 

ai = ol\, 0:3 = —4 + 2ai, 05 = 4 — 2a\ — 04, aa = —2a\ — 0:4, 
a 7 = ai + a 4 , 3 = 4 + 2 , 04 = — 4 + 0i, 

05 = As = -07 = fro = "2, 08 = 2-01, 09 = -02- (5.10) 

Substituting these solutions into the numerators, we find that the full set of Jacobi 
equations are satisfied. As a consistency check, all the numerators have the required 
power counting behavior. It also turns out that the numerators of the bubble-like 
integrals (35) — (37) are zero. 

6. There is another simple cut one can apply for the basis integrals (31) — (34), as shown in 
graph (31) in Figure 9. The cut integrand involves a three-point loop amplitude, which 
implies that the numerator should vanish under this cut. This simple constraint fixes 
two more parameters: 

1 = -ai, /3 2 = 2-^. (5.11) 

Therefore, we arrive at a solution of numerators which depends on only two parameters 
{«i,a 4 }. 

5.2 Unitarity checks 

To make sure that the solution is indeed the form factor result, we perform a large number 
of unitarity cuts. In particular a non-trivial quadruple-cut is given as in Figure 11. From the 
basis integrals, we obtain the cut integrand in terms of 192 cut diagrams from following 18 
integrals: 

{(1), (2), (3), (4), (5), (6), (7), (9), (10), (12), (13), (14), (17), (19), (21), (25), (30), (31)}. 

(5.12) 
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This should be compared with the product of trees, which involves non-trivial next-to-next-to 
MHV tree amplitudes 

r 5 r 

/ J] d\ [^ nY (-h, -h, ~h, -k, ~h) Af nV ( Pl , P2 , h,k, h, h, h) (5.13) 

i=l 

J r | JMHV ( — i 1; — \ 2% — l 3j — — / 5 ) MHV (pi,P2j ^5; ^4, ^3) ^2; ^l) 

j-N MHV^_^ ^47 MHV (p!,p 2 , ^4, h, h, h) 

jrmax-non-MHV^^ — Z 5 ) .A^ 1 ^ (px , p 2 , Z 5 , Z 4 , Z 3 , Z 2 , ^l) 

The two results are compared numerically to very high precision. We find that the remaining 
two parameters should satisfy the relation 

7 

6 = a\ + -0:4 . (5-14) 

We should emphasize that this is a very non-trivial consistency check for our result, 
considering that both the cut integrand and the product of trees are very different and 
complicated expressions. The two sides are however identical to each other if and only if 
the very simple relation (5.14) holds. This leaves one free parameter in the solution which 
is written below in full detail. The numerators of the two master integrals can be given 
explicitly as 

#13 = ~(4 ■ • Pi) - (4 • P 2 )(4 • P2) - 5(4 • Pl)(4 ■ P2) + 3(4 • P 2 )(4 ■ Pi) 

+ (pi • p 2 ) [2(4 • 4) + 2(4 • Pi) - 2(4 • P2) - (4 • Pi) + (4 • P2) + 2(pi • pa)] (5.15) 
-j(ai + 1)(4 -Pu - Pi -P2)(h ■ (7pi ~P2)) , 

N 2 i = -(4 • pi) 2 - (4 • P2) 2 - 6(4 • pi)(4 • P2) + (pi • P2) [2(4 • 4) + 4(4 • pi) + (pi • pi)\ 

+(ai + 1) [(4 • P12 - Pi ■ P2? ~ ^ (4 • (4 - P12) + Pi • p 2 )(pi • P2)] , (5.16) 

where we separate the dependence on the remaining parameter into a term proportional to 
(ai + 1) which would be zero for the possible choice of (5.18). Using the Jacobi relations 
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(e) (/) (s) W 




(0 00 (*) (0 

Figure 12. Other cuts that have been checked. A few non-planar cuts which are obtained by permuting 
the legs of tree sub-amplitudes have also been checked. 



in Appendix A, one can obtain all other numerators. They are explicitly given in Table 5, 
and are also written in the form of parameter independent terms plus a term proportional to 
(ai + 1)- Since up to now all coefficients in the numerators have been integers, it is plausible 
that 

ai + 1 = 7k keZ (5.17) 

which achieves this. Moreover, the most natural choice is a\ = — 1 which would simplify the 
result considerably. This leads to the natural value of the two parameters 

qi = -1, a 4 = 2. (5.18) 

Actually, the possibility to find only integer coefficients in the ansatz is a result in itself. 

Some other cuts have also been checked as shown in Figure 12, which are all consistent 
with the above solution. A few of them requires the same constraint (5.14), however, none of 
them provides any new constraint for the parameters. In particular, ordinary unitarity cuts 
which isolate a tree level two-point form factor cannot fix the coefficient. Therefore we finally 
arrive at a one-parameter solution of the integrand, as far as the many non-trivial constraints 
we have imposed. 
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5.3 Full result 

The explicit results including the color 10 and symmetry factors are summarized in Table 5. 
The full four-loop form factor can be obtained as 

34 
■ i 

Note that six of them have zero color factor and do not contribute to the final answer. 
However, they do play important role in the above construction 11 . 
A few comments on the result are in order. 

• All the 34 integrals have the maximal number of loop propagators. There are no multi- 
point tree legs such as the snail-like diagrams appearing in the case of four-loop four- 
point amplitudes [10]. This seems to be a special feature for the two-point form factor, 
and may be also true at higher loops. 

• Our construction is based on color-kinematic duality which does not rely on the dimen- 
sion of loop momentum variables. On the other hand, in the unitarity construction one 
usually needs to worry about cut non-constructible part such as the /j, 2 terms which 
may be missed in four-dimensional unitarity. It is reasonable to believe that our result 
should already contain such information and is true for general dimension. A full set of 
D-dimensional unitarity checks would be interesting to perform. 

• All the numerators depend at most quadraticly on loop momenta. By naive power- 
counting this implies that the four-loop result has no UV divergence below the critical 
dimension of D = 5i, as expected. Note that there could be further cancellations 
beyond this hidden in the integrals. 

• The integrand solution we have obtained still depends on one parameter. In all the 
unitarity cuts we have performed, this parameter is not fixed although as argued above 
it does have a natural value 12 . It may be uniquely fixed by some other more sophisticated 
cuts yet to be preformed, which should also provide further consistency checks for our 
result. Another important check will follow from the evaluation of the integrals and 
their infrared singularity structure [5]. 

10 To be explicit, we write the color factor for SU(N C ) gauge group. To translate to general groups one can 
rewrite our results in terms of the quadratic Casimir Ca and the quartic Casimir da as outlined in appendix 
B. 

11 As a side comment, we mention that these basis with zero color factor actually have non-trivial kinematic 
contributions in some of the non-planar unitarity cuts we have checked. Therefore they in some sense do 
contain some physical information even though they don't contribute to the final result. If an observable in 
gravity would exist, we would expect these terms to play a non-trivial role there. 

12 Although not likely, it might happen that this is a real freedom at the integrand level, which means that 
after integration the dependence will vanish. 



-27- 



Table 5. The result for the two-point four-loop form factor. 



Graph 


Numerator factor 


Color factor 


Symmetry 
factor 


(1) 


s 12 




2 


(2) 


Ni 




4 


(3) 


Ni 




4 


(4) 


Ni 


8 iV c 5 ai a 2 


4 


(5) 


A T 


a a r4 c 


8 


(6) 


-S 12 [4 • (Pl - P 2 ) + Pi ■ P2] 


4iV c 4 5 ai a 2 


2 


(7) 


-N 6 


4iV c 4 5 ai a 2 


1 


(8) 


N 6 





2 


(9) 


2(4 • Pi2)(4 • P12) - s 12 (2£ 3 ■ 4 + pi • p 2 ) 

-f («1 + 1)(4 -P12 "Pl -P2)(4 "Pl2 "Pl -P2) 


4iV c 4 5 aia2 


2 


(10) 


-Ng 


4iV 4 5 ai a 2 


1 


(11) 


N 9 





2 


(12) 


-(4 • pi) [4 • (pi + 5p 2 ) - 4 • (pi - 3p 2 )] 

+(4 • P2) [4 • (3pi - p 2 ) + 4 • (5pi + pa)] 
+(pi • pa) [24 • (Pi - P2 + 4 - 4) 
-3(4 + 4) • (Pi - P2) + 24 • (Pi - P2) + sia] 
+7(0:1 + 1)(4 • P12 - Pi • P2) 

x[4 • (7pi -p 2 ) +4 • (Pi - 7p 2 )] 


2iV 4 <5 aia2 


2 


(13) 


(4-Pi)(4-(pi + 5p 2 )) 
-(4-p 2 )((4-(3pi-p 2 )) 
-(pi - pa) [24 • (4+P1-P2) 
-4 • (pi -P2) + S12] 

-i(ai + 1)(4 -pi2 - pi -p 2 )(4 • (7pi -P2)) 


2iV c 4 <5 aia2 


1 
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Table 5 (continued). The result for the two-point four-loop form factor. 



Graph 


Numerator factor 


Color factor 


kJj 111111C LI J 

factor 


(14) 


-N 13 




1 


(15) 


-N 13 





1 


(16) 


N 13 





1 


(17) 


-(4"Pl)(4-(pi + 5p 2 )) 

+ (4-P2)(4-(3i?l -P2)) 

+( Pl • p 2 ) [2(4 • 4) + 24 • (pi - P2) 
-3£ 5 ■ ( Pl - p 2 )] +i(ai + l)x 
(4 ■ P\2 - Pi ■ p 2 )(4 ■ (7pi - P2)) 


2iV 6 4 <5 aia2 


1 


(18) 


-N 17 





2 


(19) 


(4 • Pi) [4 ■ (pi + 5p 2 ) - 4 • (pi - 3p 2 )] 
-(4 • Pa) [4 • (3pi - p 2 ) + 4 • (5pi + p 2 )] 

-(pi • P2) [24 • (pi - p 2 + 4 - 4) 
-3(4 + 4)-(pi-p 2 )] 

-\(ai + 1)(4 -pi 2 - pi -p 2 ) 

x[4-(7pi-p 2 ) + 4-(pi-7p 2 )] 




1 


(20) 







2 


(21) 


-(4 • pi) 2 - (4 • P2) 2 - 6(4 • pi) (4 • pi) 
+(pi • pi) [2(4 • 4) + 4(4 • pi) + pi • pi] 

+(ai + 1)[(4 -Pi2 - Pi -p 2 ) 2 

-f(4 • (4 -P12) +pi -p 2 )(pi -p 2 )] 


(2iV c 4 + 24iV c 2 )J aia2 


2 
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Table 5 (continued). The result for the two-point four-loop form factor. 



Graph 


Numerator factor 


Color factor 


Symmetry 
factor 


(22) 


N 21 


24 TV 2 S 


2 


(23) 


N 2 i 


24 N 2 S 


2 


(24) 


N21 


24iV 2 <5 aia2 


1 


(25) 


N 2 i 


(2iV c 4 + 24iV c 2 )5 aia2 


4 


(26) 


-(4 ■ vi? - (4 ■ P2? - 6(4 ■ pi)(4 ■ V2) 
+(4 ■ Pi) [4 ■ (pi + 5p 2 ) - 4 ■ (pi - 3pa)] 
-(4 -p 2 )[4 ■ (3pi -P2) + 4 ■ (5pi + p 2 )] 

+(pi-p 2 )[24-(pi 2 +4-4+4) 

+3(4 + 4)- (P1-P2) 

-2^4 • (p a -p 2 ) "Pl -P 2 ] 

+(ai + 1){(4-Pi2-Pi -p 2 ) 2 

-7(4 • (4 -P12) +Pi -p 2 )(pi -p 2 ) 

-1(4 -P12 - Pl -P 2 ) 

x [4 • (7pi - p 2 ) + 4 • (Pl - 7p 2 )]} 


(2iV c 4 + 24iV 2 )5 aia2 


4 


(27) 


-(4 • pi) 2 - (4 • p 2 ) 2 - 6(4 • pi)(4 • p 2 ) 
+(4-pi)(4-(pi + 5p 2 )) 
-(4-p 2 )(4 • (3pi -p 2 )) 
+(pi-p 2 )[24-(pi 2 + 4-4) 
+4 • (pi - p 2 ) - pi • p 2 ] 

+(ai + 1) [(£3 • P12 - pi ■ p 2 ? 
-7(4 -Pi2 - Pi -p 2 )(4 • (7pi -p 2 )) 

-f (4 ■ (4 -P12) +pi -p 2 )(pi -p 2 )] 


24iV 2 <5 aia2 


1 


(28) 


^27 


24 /V 2 S 


1 


(29) 


N 27 


24iV 2 <5 aia2 


1 


(30) 


N 27 


(2iV c 4 + 24iV 2 )5 aia2 


1 
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Table 5 (continued). The result for the two-point four-loop form factor. 



Graph 


Numerator factor 


Color factor 


Symmetry 
factor 


(31) 


-4(^3 -pi-pi- P2)(^3 -P2-P1- Pi) 

\ ( m 1 1 ^ J ( m m m \2 

+ {ai + i)\{l3 ■ Pu - Pi ■ P2) 

-f(4 -4) (pi - Pi) - \(h -pi2 - pi - pi) 

x [*3 I'Pl P2j "JPI P2j/ 


(2N* + 24N?)5 aia2 


1 


(32) 


^31 


24iV c 2 5 aia2 


1 


(33) 


iV 31 


24 N 2 S 


1 


(34) 


N 31 


24 AT 2 r) 


1 



6 Conclusion and outlook 



In this article the extension of ideas concerning color-kinematic duality from scattering ampli- 
tudes to form factors has been studied. The focus in this article has been the consequences of 
this extension at the loop level. In several examples it has been demonstrated that the duality 
leads to all known results in the literature up to three loops for stress-energy multiplet form 
factors in maximally super symmetric Yang-Mills theory. While these results are obtained in 
a simpler and more uniform way than the original derivations, they also provide non-trivial 
support for the conjectured duality and the validity of the construction. 

The same construction has also been applied at the four-loop level for the two-point form 
factor. A four-loop solution has been found that satisfies all color-kinematic relations and 
passes a considerable number of non-trivial unitarity checks. An interesting feature in this 
approach is that planar and non-planar integrals are obtained simultaneously. The integration 
of the resulting set of integrals is working in progress [5]. The five-loop two-point integrand 
should also be obtainable with the same methods. 

Our methods should also be applicable to form factors in less supersymmetric theories 
such as QCD. This is definitely a direction worth pursuing further as the duality relates in 
general planar and non-planar contributions. The interplay between these in the duality is 
in general a phenomenon which will be interesting to explore further. In the current form 
color-kinematic duality is only applicable to matter in the adjoint representation; a possible 
extension to matter in the fundamental representation would be very interesting as well. 

More generally, it should be fruitful to see how far and wide the general ideas of color- 
kinematic duality can be applied beyond scattering amplitude computations. These can be 
taken as a hint that there is sense to be made of the duality at a Lagrangian level in the 
gauge theory. Indeed, the existence of a Lagrangian formulation of gauge theory which yields 
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an explicitly color-dual perturbation theory would be a natural explanation of the duality 
for form factors. This has so far only been achieved for low numbers of points [12] or in 
the self-dual sector of the gauge theory [39]. Moreover, it would be interesting to find the 
observable in the gravity theory which corresponds to the 'square' of our form factor results. 
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A Four-loop Jacobi relations 

In this appendix we give a set of Jacobi relations of four-loop two-point form factor, from 
which one can generate all numerators by using the two master integral numerators: N\ 3 , N21. 

N 2 = N w [ Pl ,p 2 ,£ 3 ,£ 4 ,£ 6 , Pl2 - 4 + 4 - 4] + N w [ Pl ,p 2 ,pi2 - £3, -4+4,4,4 - 4](A.l 
JVi = N 3 = 2V 4 = N 5 = N 2 , (A.2 

N 6 = N 13 [pa , p 2 , 4, 4 , P12 - 4, 4] + A 17 [pi , P2, 4, 4 - 4 + 4, 4 , 4] , ( A.3 

iv 8 = -jv 7 = iv 6 , (A.4 
N 9 = N 13 [p 2 , pi, P12 -4,4,4,4] + iVi9bi,P2,4,4,4,W2 -4 -4] , (A.5 

iVn = -iVio = ^9 , (A.6 

n 12 = N w [p 2 ,pi,pi2 - 4,4 + 4 - 4,4,4 - 4] - iVi7bi,P2,4,4 - 4,pi 2 - 4 -4,4] : 

(A.7 

N 16 = -N u = -N 15 = N 13 , (A.8 
iVis = iV28[pi,P2,4,4,pi - 4,P2 - 4] - N2s\pi,P2,pi2 - 4,4 - 4, pi - 4,P2 - 4] , 

(A.9 

N 17 = -N 18 , (A.10 

A r 20 = N 2 8\pi,P2,Pl2 - 4, Pi - 4,4 - 4,4] - A r 28[P2,Pl,Pl2 ~ 4,P2 - 4,4,4] , (A. 11 

N 19 = N 20 , (A. 12 

N 22 = N 23 = N 2i = N 25 = N 2 i , (A. 13 

A f 26 = ^25 - A^12 (A. 14 

A f 28 = iVi 3 + N21 , (A. 15 

N 27 = N 29 = N 30 = N 28 , (A. 16 

A^si = N 13 [ Pl ,p2,£ 3 ,pi -4,4, p 2 -4] +^28^1,^2,4,4 + 4,4, p 2 -4], (A. 17 

A f 32 = ^33 = A f 34 = JV 3 1 • (A. 18 

Note that Ni = A r j[pi,p2, 4, 4, 4, 4] if not specified. 

B Classifying color factors for 2-point correlation functions to eight loops 

In this appendix the independent group theory structures are classified which appear in the 
calculation of the /-loop two-point form factor in any gauge theory with adjoint matter only. 
Maximally supersymmetric Yang-Mills theory is a particular example of such a theory. The 
main observation is that these will always be proportional to 5 ab . Contracting this with 
itself gives a vacuum graph in color space, times a universal factor which is the number of 
generators in the adjoint. In SU(N) this number is Na = (N 2 — 1). The classification in 
this appendix applies also to general two-point correlation functions for two operators in the 
adjoint representation. 
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loops 


color structures 


1 




C A 


2 




r 2 


3 




r 3 

A ~ 


4 




C\ C?44 


5 




C?44Ca 


6 




C\ d^C\ C?444 


7 


c\ 


d^C\ ^444 CU ^644 


8 




d±tJ3\ dmC\ d^CA 1 






C?844 <^4444a ^44446 A^G^ J 



Table 6. Classification of group theory invariants in a two-point calculation at a fixed number of 
loops. Tildes indicate that the 'd'-typc invariants in the output of COLOR have to be divided by a 
common factor of Na 

All vacuum graphs of a fixed number of loops consisting of only three vertices can be 
generated using DiaGen [32], excluding self energies and selecting 1PI graphs. In a non- 
supersymmetric theory the self-energy graphs are non-trivial, but will only contain color 
structures obtained at a lower loop order, multiplied with a power of Ca- 

The output of the diagram generator is parsed for each graph by Mathematica into 
a definite color trace by assigning a structure constant to each three vertex. This index 
contraction is passed to a FORM [43] program which itself uses the COLOR package [44]. 
After running the program the output of this is read back into Mathematica. This cycle is 
repeated for each graph in the output. 

The result is a list of group theory factors for each separate graph, given in table 6. 
The labels for the invariants are produced by the COLOR package. By computing a small 
choice of color traces by explicit Fierzing (which is many orders of magnitude slower than the 



COLOR package) we obtain up to seven-loop order 

C A = N (B.l) 

N A = N 2 -1 (B.2) 

J44 = — (36 A 2 + iV 4 ) (B.3) 

d~444 = ^ (324 A 2 + 135 A 4 + A 6 )) (B.4) 

d~644 = tt^t (3564 A 3 + 225 A 5 + A 7 ) (B.5) 



The form factor at I loops is proportional to g 2 m , which is usually combined with A into the 
't Hooft coupling A = Ng 2 m . In terms of A and A it is seen from table 6 that the first non- 
planar correction sets in at four loops, while the next-to-non-planar correction starts at six 
loops. This pattern is expected to continue based on the general arguments of [45]. Note that 
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these results can be used to reverse-engineer the color factor of any result given for SU (N) 
two-point functions to more general invariants. 

The number of different color structures appearing at a certain loop order for two-point 
functions comes down to classifying trivalent graphs at this loop order with two endpoints, 
modulo Jacobi relations. In particular, the two external legs can be counted as identical. 
This turns out to be a math problem which has been studied in various places, see [46] 
and references therein. This paper provides a conjecture for the number q of independent 
Casimirs appearing at loop order I as a generating function: 



Our explicit computations are in agreement with these numbers up to seven loops. The 
discrepancy at eight loops is most likely resolved by finding relations between the Casimir 
invariants COLOR introduces at this loop order: this set is not guaranteed to be minimal 
beyond seven loops by the used methods. 
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